home *** CD-ROM | disk | FTP | other *** search
/ Shareware Grab Bag / Shareware Grab Bag.iso / 004 / mathstat.arc / M19.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1982-04-14  |  1.0 KB  |  56 lines
  1. 0  CLS:PRINT
  2. 10   PRINT "REAL ROOTS OF POLYNOMIALS: NEWTON"
  3. 20   PRINT 
  4. 30   DIM A(11),B(11)
  5. 40   FOR I = 1 TO 11
  6. 50  A(I) = 0
  7. 60  B(I) = 0
  8. 70   NEXT I
  9. 80   PRINT "DEGREE OF EQUATION";
  10. 90   INPUT N
  11. 100   FOR I = 1 TO N + 1
  12. 110   PRINT "COEFFICIENT A(";I - 1;")";
  13. 120   INPUT A(I)
  14. 130   NEXT I
  15. 140   FOR I = 1 TO 10
  16. 150  B(I) = A(I + 1) * I
  17. 160   NEXT I
  18. 170   PRINT 
  19. 180   PRINT "GUESS";
  20. 190   INPUT X
  21. 200  Q = 0
  22. 210  S = 1
  23. 220  F1 = 0
  24. 230  F0 = 0
  25. 240  Q = Q + 1
  26. 250   FOR I = 1 TO N + 1
  27. 260  F0 = F0 + A(I) * S
  28. 270  F1 = F1 + B(I) * S
  29. 280  S = S * X
  30. 290   NEXT I
  31. 300   IF F1 = 0 THEN 360
  32. 310  S = X - F0 / F1
  33. 320   IF X = S THEN 380
  34. 330  X = S
  35. 340   IF Q > 100 THEN 490
  36. 350   GOTO 210
  37. 360   PRINT "DERIVATIVE = 0 AT X = ";X
  38. 370   GOTO 180
  39. 380   PRINT 
  40. 390   PRINT " ROOT"," ERROR","DERIV."
  41. 400   PRINT X,F0,F1
  42. 410   PRINT 
  43. 420   PRINT "ANOTHER VALUE (1=YES,0=NO)";
  44. 430   INPUT A
  45. 440   IF A = 1 THEN 170
  46. 450   PRINT "ANOTHER FUNCTION (1=YES,0=NO)";
  47. 460   INPUT A
  48. 470   IF A = 1 THEN 30
  49. 480   GOTO 550
  50. 490   PRINT "100 ITERATIONS COMPLETED:"
  51. 500   PRINT " X = ";X;" F(X) = ";F0
  52. 510   PRINT "  CONTINUE (1=YES,0=NO)";
  53. 520   INPUT A
  54. 530   IF A = 1 THEN 200
  55. 550   RUN "MATH"
  56.